Broken Ergodicity in classically chaotic spin systems 
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A one dimensional classically chaotic spin chain with asymmetric coupling and two different 
inter-spin interactions, nearest neighbors and all-to-all, has been considered. Depending on the 
interaction range, dynamical properties, as ergodicity and chaoticity are strongly different. Indeed, 
even in presence of chaoticity, the model displays a lack of ergodicity only in presence of all to 
all interaction and below an energy threshold, that persists in the thermodynamical limit. Energy 
threshold can be found analytically and results can be generalized for a generic XY model with 
asymmetric coupling. 

PACS numbers: 05.45Pq, 05.45Mt, 03.67,Lx 



I. INTRODUCTION 

Systems with few degrees of freedom but good chaotic 
properties, can be characterized by standard statisti- 
cal properties, for instance diffusive random-walk Uke 
behavioral or appearance of stationary distribution for 
single particle occupation numbers y, Q . 

It is also generally assumed that, due to a sufficiently 
strong interaction, chaos will provide the mechanism in 
order to have ergodicity on the energy surface. Generally 
speaking, the degree of chaos will depend on the strength 
of the interparticle interaction and on the number of par- 
ticles: typically, the larger the number of particles, the 
stronger is chaos, thus leading to more suitable condi- 
tions for ergodicity and statistical investigation. 

This common lore is far from being satisfied by most 
physical systems, as explained by R.G. Palmer in his sem- 
inal paper where, using his words: "The breakdown of 
ergodic behavior is discussed as a general phenomenon in 
condensed matter physics" . In his paper the breakdown 
of ergodicity as due to an "effective" confinement of the 
system, respect to the observational timescale has been 
put forward. Nonetheless, to the best of our knowledge, 
there are no single examples where thresholds for broken 
ergodicity can be found analytically and, more impor- 
tant, where the relation with the interaction range has 
been remarked explicitly. 

Only in recent years the relation between the dynam- 
ical properties, like chaos or ergodicity, and the kind of 
interparticle interaction has been studied |^. Indeed, 
long-range interacting systems display peculiar proper- 
ties from the point of view of statistical mechanics. It 
is well known for instance that, for such systems, the 
canonical and the microcanonical ensemble give differ- 
ent results, 0, thus questioning the validity of statistical 
mechanics. Moreover, chaos suppression in long-range in- 



teracting systems when the number of particles increases 
was recently found|3,0,Q^, and analytically investigated 
[llj l , thus leading to anomalous statistical behavior in the 
thermodynamical limit. 

In this Paper, we investigate a one-dimensional spin 
chain model with two different ranges of interaction, 
nearest neighbors and "infinite" (all particles are inter- 
acting with all the others) and we show that, even in 
presence of dynamical chaos a non-ergodic behavior 
is found, for any finite number of particles. Such non- 
ergodic behavior is mainly due to a unconnected phase 
space and the threshold for non-ergodicity can be ob- 
tained analytically. Moreover such non-ergodicity per- 
sists in the thermodynamical limit, actually the ratio be- 
tween the disconnected portion of the energy range and 
the whole one approaches one, as the number of particles 
increases. 

We have also found that, for all-to-all interaction, when 
the number of particles increases, the system dynamics 
becomes more regular. We also give numerical evidence 
that there is room for chaos in the non-ergodic region, 
for any finite number of particles. 

For sake of comparison, we also consider the same 
model with short range interaction (nearest neighbor cou- 
pling) and we prove that this kind of disconnection does 
not exist. Of course, this does not imply that the system 
is ergodic, as also supported by our numerical simula- 
tions that indicate the presence of regular dynamics in 
some region of the energy space. 



II. THE MODEL : CHAOTICITY VS 
INTERACTION RANGE 

Our model is a variant of the one-dimensional Heisen- 
berg model for iV-spins. The Hamiltonian is given by: 
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(1) 



where stands for nearest neighbor (Af- 

interaction), or infinite range couplings (aU-to-aU, A- 
interaction) and J is a positive constant. 

Hamihonian gives rise to the standard nonhncar 
equations of motion: 



(2) 



where {j) is a shorthand notation for j ^ i (in the case 
of .A-interaction) , and j = i ± 1 when A/" interaction 
(without periodic boundary conditions) is assumed. 

Constants of motion are the energy H = E, and the 
N squared moduh \si\'^ (which we set equal to 1 for sim- 
plicity). For iV = 2, an additional constant of motion, in 
involution with the others, is given by s\ — S2, and the 
system becomes exactly integrable. 

Let us notice that, in general, the only free parameters 
are the total energy E, the interaction strength J and the 
number of particles N . Moreover, for any finite number 
of particles iV, the energy is bounded \E\ < EmaxiJ, N). 
One can also give a rough estimate for such border as 
Emax ^ JN^ for yl-interaction and E„iax ^ JN for J\f- 
interaction, even if the rigorous bound will be given be- 
low. 

In order to make a fruitful comparison between Af and 
A interaction and according to a general prescription [T^ 
we rescale, for A interaction, the strength J to the num- 
ber of particles N, setting J = 2I/N. In this way the 
energy of A systems scales with the number of particles 
in the same way as for Af systems. 

For N > 2, and sufficiently strong interaction (both 
Af and A interaction) the system is chaotic in a large 
energy range \E\ < \Ech{j,N)\ < Emax, as indicated by 
a maximal positive Lyapunov exponent. From now on, 
all results will be restricted to the chaotic energy region 
only. 

Chaotic properties of the model has been anal yzed 
within the frame of Standard Lyapunov analvsis [13. 1 1 5j . 
In Fig.(^ we show the maximal Lyapunov exponent as a 
function of the energy per particle e = E/N for different 
number of particles and A (upper) and A/" (lower) inter- 
action. As one can see, while the maximal Lyapunov 
exponent Xmax depends on the energy per particle in 
some peculiar way, there is a large energy region where it 
is appreciably different from zero for both interactions. 
From the same picture it is clear that while increasing 
the number of particles in the case of A/" interaction, does 
not produces an appreciable variation of Xmax , the same 
variation in N, for A interaction, effectively decreases the 
value of Xmax ■ Even if it is difficult to give an exact scal- 
ing relation of the maximal Lyapunov exponent with all 
the system parameters {E^ iV, J) the approximate rela- 
tion Xmax[E = 0,/ = 3) - 1/iV" holds, where a = for 
N interaction and a — 1 for A interaction, see Fig. 
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FIG. 1: a) Maximal Lyapunov exponent, Xmax, for different 
number of particles A'' = 5 (black) , N = IQ (blue), = 25 
(red), N = (yellow) and A'^ = 100 (green), as a function of 
tfie energy per particle e = E/N, for ^-interaction a), and 
for A/'-interaction 6). Error bars have been obtained as the 
standard deviation from a set of 20 initial trajectories on the 
same energy surface. Trajectories have been integrated for a 
time of 10*. In a) the interaction constant has been chosen 
as J = 2I/N, with 7 = 3, while in &) J = 7 = 3. 



The numerical results and the corresponding approxi- 
mate relations are shown in Fig.© and indicate the ab- 
sence of a chaotic dynamics for all-to-all interacting sys- 
tems with a large number of particles. 

The general trend where the maximal Lyapunov expo- 
nent decreases with the number of particles can be un- 
derstood from the fact that in the limit N ^ 00 and for 
A interaction the model becomes close to an integrable 
one. 

Indeed the Hamiltonian per particle can be written as: 



H I 

= Ihn H h^ 

N ' N 



(3) 



where ho 



2-m2and/.i = iEti [(^r)'-«n 

where we have defined the average magnetization mj, = 

While both ho and hi remain on order 1 in the limit 
of large N, and fixed interaction strength I, the constant 
in front of hi goes to zero and the /iq term dominates. 
A close inspection [l6j. and our numerical results indicate 
that ho represents an integrable model for any choice of 
parameters. 




FIG. 2: Maximal Lyapunov exponent as a function of the 
number of particles, for the same parameters of Fig.Q and 
E = 0. Triangles refer to A interaction, while crosses are 
for J\f interaction. Also shown the best fit lines Xm = 1.1 ± 
0.2(blue) and \a = (4.1 ± 0.1)/iV (red). 



III. ERGODICITY : NUMERICAL RESULTS 

While in two degrees of freedom systems the absence of 
ergodic motion in presence of chaos is somehow obvious 
(due to the presence of invariant tori), for many degrees 
of freedom systems the same occurrence is far from being 
trivial. In spite of that, we will now show that, for A 
interaction, an energy threshold exists below which one 
cannot have ergodic motion. 

Let us first evaluate ensemble and time averages for the 
y magnetization and from them, respectively, the proba- 
bility distributions PP{niy) and P*{my). 

Phase and time distributions of the mean magnetiza- 
tion are shown in Fig.Q for the A/'-interaction (right 
column) and y^-interaction (left column) with the same 
number of particles. All cases are characterized by a 
chaotic dynamics, as given by a maximal positive Lya- 
punov exponent, see Fig. 

As one can see, while for nearest neighbor interaction 
there is a good correspondence between the two aver- 
ages (this of course does not mean ergodicity), in the 
case of all-to-all interaction there are strong deviations, 
in the lower energy case, see Fig.(|2h)- More precisely 
one trajectory with an initial ruy > 0(< 0) cannot reach 
a region with my < 0(> 0), below some energy thresh- 
old. Another difference is related to the presence of a 
two-peaks distribution: this indicates the presence of a 
phase transition and it will be the subject of a separate 
paper. 

Strictly speaking, from these numerical results one can 
only infer that the time of transition from one peak to 
the other one is much larger than the simulation time. 
Nonetheless, in the next section we will show that such 
transition time does not exist. The origin of this devia- 
tion is that for those energy values less than the energy 
threshold the phase space is metrically decomposableflTj 
thus making the A system rigorously non-ergodic. 
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FIG. 3: Distribution of mean magnetization along the y- 
axis for A interaction (left column) and Af interaction (right 
column). As black lines we show the phase distributions, 
P^{my), while red ones stand for time distributions, P^iniy). 
Phase distributions have been obtained taking M — 10^ dif- 
ferent initial conditions upon the energy surface E — —5, a) 
and b) and E — —1, c) and d), while time averages have been 
obtained by computing one single trajectory in time on the 
same energy surface for a time T = 10^ and taking values at 
time steps At = 1. Results have been checked on increasing 
the number of initial points upon the energy surface for the 
phase average, and on increasing the integration time T for 
the time average. Here is A'^ = 10, and / = 3. Maximal Lya- 
punov exponents are positive for all the cases shown in this 
figure, as can be deduced from Fig.Q. 



IV. NON-ERGODICITY BORDER : 
ANALYTICAL RESULTS 

The existence of a threshold for the non ergodicity in 
case of A interaction can be proved by the following argu- 
ment. Let us assume that initially the spins are oriented 
in such a way that ^^(0) > with the energy E. It is 
clear that if we are able to prove that when E < Ene one 
cannot have the solution niy = on the energy surface 
H = E, it means that the trajectory should keep, for any 
time my{t) > 0. Since reversal of y single spin compo- 
nents results in change of y magnetization but not of the 
energy, the same argument can be applied to my{0) < 
as well. In other words the space is metrically decom- 
posable, since one trajectory with my > cannot reach 
the region where my < (and viceversa), although both 
regions belong to the same energy surface. 

In order to find the explicit value of Ene one should 
consider the problem of minimizing the value of the en- 
ergy E under the constraint my — {l/N) X^iLi — 0- 
principle this task can be faced by means of the Lagrange 
multipliers. Nevertheless from a very simple argument we 
can estimate the critical energy Ene- 

Let us first derive the minimal energy, Emm , and then 
the non-ergodic energy, Ene, for both Af and A systems. 
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A. Minimal energy 



B. Non— ergodic threshold 



We can find the minimal energy by minimizing the 
Hamiltonian ^ under the constraint: 



l,...,iV (4) 



We can take conditions Q into account, switching to 
spherical coordinates: 



sin(6'i)cos((/>i) 
sin(6'i)sin(0j) 
cos(6'i). 



(5) 



We can now rewrite Hamiltonian as: 
H = J ^ shi{0i) am{Oj) cos((/)i + 

where 



(6) 



J = 



(2J/7V) for ^-interaction 
/ for TV— interaction 



(7) 



One absolute minimum value is obtained, for instance, 
setting: 



sin(0,) = 1 

COs((/)i + = —1 



(8) 



where j = i + 1 for Af system and j ^ i for A system. 
While the first equation of ||SJ) can be easily satisfied 
by both Af and A system, assuming 



TT 

2 



(9) 



some care should be taken in order to fulfill the second 
one. 

For TV system, we can put for instance: 



for i even 
TT for i odd 



(10) 



On the contrary, for A system, we set 



= 7r/2 for i = l,...,N. 



(11) 



Taking into account that there are — 1 couplings 
in the TV system and N{N — l)/2 in the A system, the 
minimal energy is given by (for both interactions): 



Ermn^ -HN -I). 



(12) 



We define the non-ergodic threshold Ene in such a 
way that for any energy E < Ene the system is met- 
rically decomposable. In order to find this non-ergodic 
(or disconnection) energy we have to find the minimum 
of Hamiltonian under the additional constraint 



rriy = 0. 



(13) 



For A system Hamiltonian (^) can be written as, 



H 



TV 



''^^f + j;^T.(''ir + INml-INml (14) 



Let us now search for the minimum of H14(l under the 
constraints Q and (|13l) . In 114|) the only term that can 
be negative under the given constraints is the first one, 
therefore. 



E„ 



> 



If it is possible to minimize this term putting at the same 
time all the other to zero, we get the non-ergodic energy. 

When TV is even we can put half spins with sf = -1-1, 
and half such that sf = — 1. In this way nix = 0, sj' = 
and niy = 0, so that: 



E„ 



L for TV even 



When TV is odd we cannot minimize the first term in 
(|14|l . and at the same time let the other terms to be zero, 
so that: Ene > —I , for TV odd. 

Anyway it is easy to give an upper bound for Ene which 
is sufficient for our scope. Indeed we can arrange TV — 1 
spins as in the previous case (TV even), and for the last 
spin we assume = ±1. In this way we have rUx = 0, 

= and ruy = 0, so that: 

-/ < Ene <-! + I/N, for TV odd 

For TV system, the non-ergodic energy, if defined as 
the minimum energy for which Ea. H13|l holds, coincides 
with the minimal energy. Indeed the same values of (jji 
(|10|l and 9i which minimize the Hamiltonian, satisfy 
condition H13|l too, so that: 

Ene = Emin = —/(TV — 1) 

This of course does not exclude the possibility to have 
some other constraint that produces a disconnected phase 
space (even if such conclusion is not supported by our 
numerical simulations). 

In Tables III and Hll we summarize the minimal and non- 
ergodic energies found. 

As an interesting consequence of our results we note 
that the ratio: 



\Ene En 
lEm 7.r). 



(15) 
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-I{N - 1) 


-I 


Nodd ■■ 


-I{N-1) 


(-/, -/ + I/N] 



TABLE I: Minimal and non-ergodic energy for A system. 
Note that in the case of odd number of particles we give a 
lower and an upper bound for the non-ergodic energy. 



From this equation we can estimate, for large N 



J , 



As for the non-ergodic energy, following the same proce- 
dure described in the previous section, we have: 





Emin 


Ene 




-I{N-l) 


-I{N-1) 



TABLE IL Minimal and non-ergodic energy for M system. 



En 



so that 



-'^N'^rj for rj>Q 
for 77 = 

^Nrj for ?7 < 



(18) 



between the disconnected portion of the energy range and 
the total energy range, in the limit of a large number of 
particles at fixed interaction strength, goes to one for A 
interaction, and is exactly zero for M interaction, thus 
showing that in this limit the energy range is completely 
disconnected for A interaction only (note that in the pos- 
itive region of the energy range the same argument used 
to show the existence of a non-ergodic energy can be ap- 
plied considering nix instead of my). 

As one can see, for A interaction, as the number of 
particles increases the region of non-ergodicity does not 
decrease thus indicating that caution should be taken 
when ergodicity is tacitly assumed due to small interac- 
tions between many particles. 

The presence of the non ergodic region is not related 
with the chaoticity of the system. Indeed, since the Lya- 
punov exponent increases linearly with the interaction 
strength, /, for any N we can find an If. value such that 
there is well developed chaos for E < Ene and I > Ic, so 
that the system is chaotic and non-ergodic, below Ene- 



GENERALIZATION 



1 — r/ for 77 > 
1 for r/ < 



(19) 



As one can see r — only for r/ = 1, while r 7^ for 
asymmetric coupling. 

For the M case, we have that 



En 



~J{N-1), 



while 



~JiN - 1) < Ene < -J{N ~ 3) 

as one can see setting, for instance, 

J. J 1 for 7 odd 
' 1—1 for i even 



(20) 



Note that for iV ;» 1, Emin — Ene so that r = for 
any rj. 

Thus, we can conclude that broken ergodicity is generic 
for asymmetric coupling and all-to-all interaction. 



Let us finally mention that the previous results can be 
easily generalized. 

Let us consider a generic XY Model: 

H=-jJ2{s!sy + rjs^s^) (16) 
{hi) 

where jryj < 1 and J > 0. Following the same arguments 
given above, we can show the existence of a non-ergodic 
energy and roughly estimate their dependence on system 
parameters. 

We will mainly focus on the ratio between the discon- 
nected portion of the energy range and the total energy 
range, r, Eq. (|15|l . in the limit of a large number of par- 
ticles. 

For A interaction, we can rewrite the Hamiltonian Hl()|) 
as: 

^ = ^ E [^(O' + (-n^] - {vml + ml) (17) 



VI. CONCLUSIONS 

In conclusion we have studied a spin chain model from 
the dynamical point of view for two different ranges of 
inter-spin interaction (nearest neighbor and all-to-all). 
Both models are mainly chaotic, in the sense that the dy- 
namics is characterized by a positive maximal Lyapunov 
exponent in a wide energy region. Nevertheless, in the 
case of an infinite range of interaction, the motion does 
not explore the whole energy surface. This can be un- 
derstood on the basis of a special topology of the phase 
space. An analytical estimate gives an energy thresh- 
old below which the system is rigorously non-ergodic: 
the existence of the threshold has been also confirmed 
by numerical simulations. Moreover, for any interaction 
strength, on increasing the number of particles the range 
of energy values characterized by such non ergodic mo- 
tion spans a finite portion of the whole energy range, if 
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the interaction has an infinite range. This could be par- 
ticularly relevant in the study of multidimensional sys- 
tems, where usually ergodicity is tacitly assumed for a 
large number of weak-interacting particles. We have also 
shown that this phenomenon is generic for XY model 
with asymmetric coupling. 

Further studies about the consequences of broken er- 
godicity with respect to phase transitions, and the com- 
parison with the correspondent quantum system are un- 



der current investigations'!' 
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